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ABSTRACT 

Context. The Serpens North Cluster is a nearby low mass star forming region which is part of the Gould Belt. It 
contains a range of young stars thought to correspond to two different bursts of star formation and provides the 
opportunity to study different stages of cluster formation. 

Aims. This work aims to study the molecular gas in the Serpens North Cluster to probe the origin of the most 
recent burst of star formation in Serpens. 

Methods. Transitions of the C"0 and C'^0 observed with the IRAM SOm telescope and JCMT are used to study 
the mass and velocity structure of the region while the physical properties of the gas are derived using LTE and 
non-LTE analyses of the three lowest transitions of C'*0. 

Results. The molecular emission traces the two centres of star formation which are seen in submiUimetre dust 
continuum emission. In the ~40 Mg NW sub-cluster the gas and dust emission trace the same structures although 
there is evidence of some depletion of the gas phase C"*0. The gas has a very uniform temperature (~10 K) 
and velocity (~8.5 kms"') throughout the region. This is in marked contrast to the SE sub-cluster. In this region 
the dust and the gas trace different features, with the temperature peaking between the submiUimetre continuum 
sources, reaching up to ~14 K. The gas in this region has double peaked line profiles which reveal the presence 
of a second cloud in the line of sight. The submiUimetre dust continuum sources predominantly appear located in 
the interface region between the two clouds. 

Conclusions. Even though they are at a similar stage of evolution, the two Serpens sub-clusters have very different 
characteristics. We propose that these differences are Unked to the initial trigger of the collapse in the regions and 
suggest that a cloud-cloud collision could explain the observed properties. 
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1. Introduction 

Despite the importance of understanding the processes 
driving the formation of stars in the Galaxy, little is 
known about the role played by molecular cloud kine- 
matics on triggering or suppressing star formation. 
Since most stars form in clusters (Lada & Lada 2003), 
the kinematics of young stellar clusters in which the ini- 
tial conditions of clustered star formation are still im- 
printed in the gas and dust emission properties can pro- 
vide important insights into the dominant mode of star 
formation (e.g. Peretto et al. 2006). 
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One such young and nearby cluster is in the Serpens 
Molecular Cloud (MC). Located at -260 pc (Straizys 
et al. 1996), the optical extinction map of the cloud 
covers more than 10 deg^ (Cambresy 1999). However, 
the majority of the star formation is occuring in three 
clusters covering approximately 1 .5 deg^ (Enoch et al. 
2007). The most active region is the Serpens Main 
Cluster (hereafter Serpens) which has a surface den- 
sity of YSOs of 222 pc~^, compared to 10. 1 pc"^ in the 
rest of the Serpens cloud (Harvey et al. 2007a). In this 
main cluster, the average gas density is around 10* cm~^ 
(Enoch et al. 2007) with H2 column densities greater 
than 10^^ cm"^ in the cores. The high density of proto- 
stars in this main cluster seems to indicate an early stage 
of evolution where the cluster gas may still be infaUing 
into the cores (Williams & Myers 1999; 2000;Hurt et al. 



Table 1. Submillimetre sources in Serpens Main Cluster 
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1996). The star formation rate in this main cluster is 
56 MoMyr 'pc^^, ~ 20 times higher than in the rest of 
the cloud (Harvey et al. 2007a). 



Amongst the youngest YSOs found in Serpens there 
are ten Class and I protostars which are detected in 
850 nm dust continuum emission (e.g. Hurt & Barsony 
1996, Davis et al. 1999), hereafter referred to as submil- 
Umetre sources (shown on Fig. 1, Table 1 and discussed 
on Sec. 2.2). These are distributed within ~ 0.2 pc^ 
and divided between two sub-clusters, one to the north- 
west (NW) and one to the southeast (SE). These sub- 
millimetre sources power a number of outflows, which 
have been studied using several difl'erent approaches 
(e.g. Eiroa et al. 1992; Davis et al. 1999; Hodapp 1999; 
Davis et al. 2000; Graves et al. 2010). Figure 1 also 
shows the Spitzer MIPS 24 //m emission tracing the 
young protostars classified as Class I or 0. The old- 
est objects in the area shown on the image are a few 
flat spectrum sources (Harvey et al. 2007a; Kaas et al. 
2004). The presence of Class II and Class III objects 
(not shown in Figure 1) dispersed over a larger area sug- 
gests that the region has undergone two episodes of star 
formation. The first, responsible for these dispersed pre- 
main sequence stars (the Class II and III sources), oc- 
curred about 2 Myr before the most recent burst which 
formed the submillimetre and 24/im protostars (Class 0, 
I and flat spectrum sources), 10^ yr ago (Harvey et al. 
2007a; Kaas et al. 2004). 



This paper focuses on the dynamical and physi- 
cal properties of the gas in Serpens using CO isotopo- 
logues observed with the IRAM 30m telescope and 
with JCMT, to probe the current properties of the re- 
gion as well as investigate any link back to the initial 
conditions under which the most recent burst of star 
formation in Serpens took place. Section 2 presents the 
observations, as well as the data reduction and analysis 
methods and techniques used in this study. In Section 
3 the structure of the gas is discussed while Section 4 
discusses its dynamics. In Section 5 its physical proper- 
ties are analysed. These results are drawn together and 
a scenario for the origin of the star formation in Serpens 
described in Section 6. 




RA (J2000) 



Fig. 1. Map of the SCUBA 850 yum continuum emis- 
sion in contours showing the position of the submillimetre 
sources (labeled). Contours at 0.4, 0.6, 1, 1.4, 1.8, 2.4 and 
5 Jy beam"' . In grey scale are the Spitzer MIPS 24 pm sources 
(Harvey et al. 2007b). All the sources seen on this figure are 
classified as being young protostars, mostly Class O and Class 
I sources, with a small number of flat spectrum sources. 

2. Data and Analysis Techniques 

2.1. IRAM Observations 

The Serpens region was observed in the J=1^0 and 
J=2— »1 transitions of C'**0 and the J=l— >0 transition 
of C"0 with the IRAM 30m telescope, using the fa- 
cility receivers, in May 2001. The observations con- 
sisted of on-the-fly maps of the region, centered at RA 
= 18''29"'57.91' and Dec = 1°12'25.2" over an area of 
approximately 10.5 arcmin-^, ~ 3' in Right Ascension 
and 3.5' in DecUnation. 

The C"0 J=1^0 data, observed at 112.359 GHz, 
have a spatial resolution of 22", a velocity resolution 
of ~ 0.052 kms"' and a noise level of ~ 0.45 K (in 
T^) in the raw map - low enough to allow the detec- 
tion and identification of the hyperfine components of 
the J=1^0 fi-ansition of C"0. C^^O was observed wifli 
spectral resolution of ~ 0.053 kms"' at 109.782 GHz 
and 219.816 GHz and with spatial resolution of 22" and 
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11" for the J=1^0 and J=2^1 transition, respectively. 
Both emission lines are detected with a good signal to 
noise, both with a one sigma noise level of ~ 0.45 K in 

The beam and forward efficiencies of the IRAM 
30m telescope (Bgff and Feff respectively) are given 
on the telescope website. From these we estimate for 
both and C"*0, a Feff=0.95 and Beff=0.72, for the 
J=l— >0 transition, and a Feff=0.91 and Beff=0.54, for 
the J=2— >1 transition. 

The main data reduction was performed using 
GILDAS software (CLASS90 and GREG). This in- 
cluded the baseline corrections, hyperfine/gaussian fit- 
ting of the data, and construction of the datacubes. 
Given the good quality of the data, the baselines were 
well fitted by a simple first degree polynomial function. 

2.1.1. 

The C'^O J=l— >0 line comprises three, partially 
blended, hyperfine features. By fitting the hyperfine 
structure (HFS) of the spectrum, the line width, velocity 
and optical depth (r) can be extracted. The line shape in 
the presence of hyperfine structure can be described by 

T(v)^Ts(l-e-'^% (1) 



where 




r(v) is the line brightness temperature, Ts is the source 
temperature and r is the optical depth. The optical depth 
is the sum over the three hyperfine components of the 
transition with r, and Vo,, - the relative weight and the 
central velocity for each hyperfine component - being v 
the velocity, cr the velocity dispersion and tq the total 
optical depth common to the three components (Fuller 
& Myers 1993). The spacing and weight of the hyper- 
fine components were adopted from Ladd et al. (1998). 
Further details about the hyperfine structure fitting pro- 
cedure in the GILDAS software can be found on the 
IRAM website'. 

To fit this hyperfine structure, the individual spec- 
trum at each pixel in the image was extracted from the 
datacube and was fitted using the procedure described 
above. A model Gaussian spectrum for each pixel was 
then reconstructed using the derived values (the peak 
intensity, line width and central velocity). Only pixels 
where both the line width and line peak intensity were 
determined with a signal to noise ratio of 5 or greater 
were considered. 

The initial fitting showed that within the uncertain- 
ties, all the emission was consistent with being optically 
thin. Therefore, to reduce the uncertainties on the fit- 
ted quantities, the HFS fitting was redone fixing the r 
at 0. 1 for the whole map, consistent with optically thin 
emission. In the final C'^O J=l— >0 modelled datacube 

' http://www. iram.fr/IRAMFR/GILDAS/doc/html/class- 
html/nodeS. html 



(Fig. 2) we were able to identify clear peaks at diff'erent 
velocities and positions in the region. A detailed study 
of these peaks is presented in Section 3.1. 

2.1.2. C'^O 

The Serpens emission is not affected by the out- 
flows in the region (Section 4. 1) and the C'^O lines have 
no hyperfine structure. Therefore, on some regions of 
the mapped cloud such as the NW sub-cluster (Fig. 1), 
lines are well represented by a single Gaussian (Fig. 3 
left panel). However, in the SE sub-cluster, the line pro- 
file has two clear peaks (Fig. 3 lower panel). The low 
optical depth of the C'**0 emission (Section 2.1.1 and 
5) together with double peaked lines in this region in 
other optically thin tracers such as N2H^ J=l— >0 (Olmi 
& Testi 2002) suggests that these two peaks trace two 
clouds along the line of sight towards this sub-cluster 
In the weaker C'^O emission with its blended hyperfine 
structure, these two velocity components are too diffi- 
cult to separate. However, since it is optically thin, the 
C'^O emission is still a reliable tracer of the total col- 
umn density. 




Fig. 3. Observed C'**0 J=l->0 spectra (black solid line), 
smoothed in velocity to 0.1 kms"' width channels, at two dif- 
ferent positions: on SMMl (top) and on SMMU (bottom); 
with the respective Gaussian fit (red line). These are examples 
of a single peaked spectrum as seen in the NW sub-cluster and 
a double peaked profile, as seen in the SE sub-cluster. 



2.2. JCMTdata 

Our excitation analysis (Section 4) makes use of JCMT 
HARP data from the Gould Belt Survey (GBS) at 
JCMT (Graves et al. 2010; Ward-Thompson et al. 
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Fig. 2. Contour maps of C'^O J=l— >0 modelled emission from the hyperfine fit, overplotted on SCUBA 850yum emission in 
gray scale. The submillimetre sources are identified on the left figure with labels and triangles. These symbols will be used 
from this point forward. Each figure represents a mean intensity map over 1 kms"' intervals, from 6 to 10 kms"'. Contours 
range from 0.4 Kkms"' increasing by steps of 0.1 Kkms"' (in T^ scale). The figure shows the SE sub-cluster with three main 
peaks: one above SMMU, one between SMM2 and SMM6, and finally, one west of SMM4. While the NW sub-cluster also 
contains three major peaks: one on SMMl, one close to SMM9, and another one south-west of SMMl. 



2007). The data used is C'**© J=3^2, at 329.330 GHz, 
with 0.055 kms"' spectral resolution, and 14" spatial 
resolution. The telescope main beam efficiency at this 
frequency is 77n,b = 0.66 (Curtis et al. 2009), and the 
rms level achieved is of the order of 0.2 K (Ta). A full 
description of these data is given in the GBS Serpens 
First Look paper (Graves et al. 2010, hereafter referred 
to as SFLPaper). 

The submillimetre continuum data at 850 fim was 
observed with SCUBA at the JCMT, with a beam size 
of 14". The initial reduction, analysis and discussion of 
these data was presented by Davis et al. (1999), where 
they estimate the overall dust properties and charac- 
teristics of the cloud (see Section 1). We have used 
the pipeline reduced SCUBA data from the Canadian 
Astronomy Data Centre (CADC) archives^ to investi- 
gate the structure of the dust continuum emission and 
for comparison with the IRAM 30m C'^0 and C'^O 
data. 

An initial inspection of the SCUBA data indi- 
cated good agreement in the source positions for those 
sources where Davis et al. determined positions from 
this same SCUBA data (SMM8 and SMMll) and as 
well as for SMM3. However, in agreement with inter- 
ferometric continuum observations (Hogerheijde et al. 
1999), the positions of some of the remaining SMM 
sources needed to be revised compared to those listed in 
Davis et al. (1999) with absolute offsets from the pub- 
lished positions greater than 5" for SMM2 and SMM6. 
Table 1 presents redetermined positions for all sources, 
extracted from the 850 fim map of Serpens, which now 
agree within 1" of the positions in the SCUBA cores 
catalogue published by Di Francesco et al. (2008). We 
estimate that these positions are accurate within the 2" 
SCUBA pointing eiTors (Davis et al. 1999). The offsets 
in RA and Dec between the revised positions and those 

- http://www.cadc.hia.nrc.gc.ca/jcmt/ 



previously published (listed in Davis et al. 1999) are 
also shown on Table 1 . 



3. Gas structure of the cloud 

To determine the structure of the molecular gas we have 
carried out a clumping analysis in 2D and 3D. Using the 
velocity information from the gas emission it is possi- 
ble to identify the individual clumps within the cloud. 
These are compared to the structure visible in the dust 
continuum. We use this analysis to quantify the sizes 
and masses of molecular gas associated with protostars, 
and carry out a virial analysis to determine the clump 
stability. 

3.1. C^''02D-Clumps 

Initially we manually extracted the small scale molec- 
ular structures for comparison with the dust seen in the 
SCUBA map. This was done based on a visual inspec- 
tion of both channel and integrated intensity maps (e.g. 
Fig. 2). The C'^0 J=l— >0 channel maps show a signifi- 
cant number of emission features which are not directly 
associated with the SCUBA cores. For this reason we 
call these molecular structures "clumps" although this 
term is often used to describe parsec-scale structures 
(Blitz 1993). 

The properties of each identified clump (Fig. 4) 
was subsequently extracted using the IDL 2D version 
of the source extraction clumpfind algorithm code by 
Williams et al. (1994) on maps integrated over the ve- 
locity range in which each clump appeared. 

With the size and the integrated intensity corrected 
for telescope efficiency, we estimated the column den- 
sity and mass of each clump (Met) assuming a tem- 
perature of 10 K, a mean molecular weight of 2.33 
and a C'^O fractional abundance with respect to H2 of 
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Table 2. Properties of the 2D-clumps 
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Fig. 4. SCUBA map of the 850 yccs\ continuum emission 
(colour scale) overplotted with the positions of the submil- 
limetre sources (triangles), the C'^O J=l-0 2D-clumps (solid 
white contours and letters) and the NW and SE sub-clusters 
(black dashed contours). The solid white contours are the in- 
tensity integrated over the correspondent velocity range for 
each clump. These are stepped by 0.2 Kkms"', except for 
weaker 2D-clumps G and D, stepped by 0.1 Kkms"' (in T^. 
The lower level of each clump is the same as the specified on 
Table 2. The NW and SE sub-clusters contours are integrated 
intensity over the entire velocity range. Contours are stepped 
by 0.25 Kkms"', with the lower contour at 1.25 Kkms"' (in 

4.7 X 10"** (Frerking et al. 1982; J0rgensen et al. 2002). 
We also calculated the clumps virial masses using Eq. 3, 
where Myh- is the virial mass, 0-^^% is the observed ve- 
locity dispersion, G is the gravitational constant and a 
is a coefficient function of the adopted density profile: a 
is 3/5 for a uniform density, 2/3 for a profile as p oc r"', 
3/4 when p oc r^-^ , and 1 when p oc f^. 



The listed virial masses of the clumps adopt a den- 
sity profile of p oc r"^. The velocity FWHM of each 



clump was estimated by averaging all the spectra as- 
signed to that clump and has an estimated uncertainty 
of ~0. 1 kms The clumps identified by this method 
will be referred to as the 2D-clumps hereafter. 

The virial mass (MvMai) and the gas mass (Mcf) 
were also calculated for the two sub-clusters, NW and 
SE. The method was the same as for the clumps except 
the density profile for the sub-cluster gas was assumed 
to be p oc r"'^, which is expected to be more appro- 
priate for these larger size regions. If the same p oc 
as for the clumps had been adopted, the derived sub- 
cluster mass would be a factor of 25% smaller 

The clumps are shown on Fig.4, and the physi- 
cal parameters summarized on Table 2, where: RApeak 
and DeCpeak are the position where the emission peaks 
within each clump; Vpeak is the velocity at the peak po- 
sition, with an uncertainty of 0.05 kms"'; area is the 
surface in the map occupied by each clump; Mcf is the 
mass of the clump calculated from clumpfind outputs; 
Mviriai is the virial mass; Ratio Mvidai/Mcf is a measure- 
ment of how bound each clump is - a gravitationally 
bound structure should have a ratio around unity, but 
given the uncertainties of these calculations, we con- 
sider a structure to be unbound if the ratio is above 
2; Iiow represents the lower contouring level assumed 
when running the algorithm for each different clump 
(increasing with steps of 0.10 Kkms"'); and, finally, 
Ipeak shows the integrated intensity in T^as measured at 
the peak position. 

Observational sources of uncertainty include the 
distance to Serpens and the line width. Uncertainties 
on the line width in particular might be a special issue 
in the SE region where the two velocity components 
observed in C"*0 may become important in broaden- 
ing the line. Systematic uncertainties in Mcf in- 
clude uncertainty in the adopted gas temperature and 
fractional abundance of C'^O. Finally, the systematic 
uncertainties on the Mviriai include source geometry ef- 
fects and the neglection of additional terms in the virial 
equation (due to external pressure, magnetic pressure, 
etc.). Amongst all the possible sources of uncertainty, 
the greatest is likely to be the factional abundance of 
C'^O, given that our non-LTE study of C"*0 at 8 po- 
sitions (Sec. 5) show a mean depletion factor of 2.5 
(App. B). Given the observational and possible sys- 
tematic uncertainties on the calculations, the virial ra- 
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tio is perhaps best seen as a useful tool to compare the 
different structures within a cloud rather than absolute 
measure of the gravitational equilibrium of any given 
clump. 

The NW and SE sub-clusters are extended regions, 
and therefore, the peak positions and velocities cor- 
respond to one of the smaller identified clumps lying 
within the sub-cluster. The NW sub-cluster peaks at the 
position of clump A (& SMMl) and the SE sub-cluster 
peaks at the position of clump F (north of SMMll). 
Similarly, the velocities quoted for the peak for the sub- 
clusters are not the mean velocity of the sub-clusters, 
but the velocity at the peak of the strongest clump. 

Note that even though both sub-clusters, SE and 
NW, have similar masses (Mcf), they each have a dif- 
ferent equilibrium status, with a factor of 3 difi"erence 
between their respective virial ratio. Interestingly, even 
when considering some depletion (Appendix B), the 
SE region is likely super-virial whereas the NW, due 
to its smaller line width, is marginally sub-virial. About 
67% of the mass in the NW region and 40% of the mass 
in the SE region is associated with the clumps. Four of 
the clumps (A, B, C and F) are individually within a fac- 
tor of three of being in virial equilibrium. Accounting 
for some depletion of C^**0 (Appendix B), Mcf could 
increase up to a factor of 2.5, which would make all 
these four clumps relatively bound structures. Clump D 
is a factor of ~ 6 super-virial, and is Ukely to be less af- 
fected by depletion as the dust densities are lower, likely 
identifying this clump as part of a more diffuse region 
which is less bound than the NW sub-cluster. Finally, 
even accounting for possible depletion, clumps, E and 
G, with mass ratios of ~ 13 and ~ 5, are likely unbound 
structures. They may either represent shocked regions 
where the line width is intrinsically high (from 1.6 to 
2.2 kms"') or regions where, as mentioned in section 
2.1.2, two blended velocity components contribute to 
the emission along the same fine of sight. The HFS fit- 
ting of a single component in this case would result in a 
broadening of the line width due to blending of the two 
components. However, in Sec. 4.3 we see that even af- 
ter separating the components, the line width of at least 
one of them is still broader than seen anywhere in the 
NW, pointing to a genuine broad line emission. In sum- 
mary, the SE sub-cluster is much more dynamic than 
the NW, with a kinematic support a few times higher, 
both when comparing individual clumps and the over- 
all sub-clusters. 

3.2. C^'' O 3D-Clumps 

Although the 2D clumpfinding is valuable for compar- 
ison with the dust continuum, it is limited in its abil- 
ity to represent the true structure of the cloud. The 
3D CLUMPFiND automatically studies the datacube in all 
three dimensions of space-space-velocity. In particular, 
3D CLUMPFIND should provide a better understanding of 
the cloud's structure where clumps may overlap along 
a line of sight but have different velocities, or where 
the emission is narrow in velocity making it weak in 



integrated intensity maps. Therefore we have comple- 
mented the 2D study of the structure of the C'^O us- 
ing the 3D version of the clumpfind within the Starlink 
package. 

Similarly to Pineda et al. (2009), we also found that 
the results on the 3D clumpfind analysis to be very sen- 
sitive to the parameters used, especially in characteris- 
ing the weaker emitting regions. Stronger clumps were 
unequivocally detected with a wide range of parame- 
ters, but changing the step size and/or the number of 
pixels per clump allowed to be adjacent to a bad pixel 
would result in the merging of several clumps into one, 
or unrealistic extensive splitting of clumps into several 
small structures, or even non-detection of some struc- 
tures expected to be detected. For this reason, the initial 
2D study is essential as a reference point to understand 
the main structure of the cloud, which could be signif- 
icantly misrepresented by relying, uncritically and ex- 
clusively on the 3D clumpfind analysis. The best config- 
uration parameters we found for this analysis were: the 
first contour level, Tiow, of 0.6 K; the global noise level 
of the data, r.m.s., of 0.2 K; and the spacing between 
the contour levels, AT, of 0.05 K. 

This analysis identified a total of 16 clumps which 
will be called the 3D-clumps hereafter. These clumps 
are shown on Figure 5 as integrated intensity maps 
in T^.kms"' plotted over the continuum 850 /im data 
from SCUBA. Table 3.2 shows the properties of the 3D 
clumps as numbered and plotted on Figure 5. The nine 
first columns are as in Table 2 the last column being the 
intensity at the peak position in T^. Once again, masses 
were calculated after correcting for the IRAM 30m tele- 
scope efficiency for the C'^O J=l— >0. 

Due to the difficulty in interpreting partial spectra 
split by 3D clumpfind between multiple spatially coin- 
cident clumps, the mean line width of the 3D-clumps 
was recovered using a different approach to that used 
for the 2D-cIumps. The velocity dispersion, cr, of each 
clump was estimated by determining the velocity range 
where the emission of the clump was above e"'''^ of 
its peak intensity. This was done by visually inspect- 
ing these thresholded channel maps of each clump. The 
quoted FWHM is 2.35 cr and has an estimated uncer- 
tainty of 0.1 kms \ twice the uncertainty of the peak 
velocity, 0.05 kms"^ 

The mass of the clumps within the NW sub-cluster 
inferred from the size and integrated intensity of 3D- 
clumps 1, 2, 5 and 8 correspond to about 65% of the 
total mass of that sub-cluster. Including clumps number 
11, 13 and 15 in this calculation, the fraction of gas- 
mass in the clumps rises to 70%. The 3D-clumps 3, 6, 
10, 12 and 14 constitute 40% of the mass of the SE 
sub-cluster. This is consistent with the results from the 
2D-clumps. 

Since the velocity structure in the region can af- 
fect the deduced clump structure, we also experimented 
with 3D clumpfind on the C'^O J=I^O data. The re- 
sults from this differed from those of C'^O only in 
that two clumps (3D-clumps 3 and 10) were subdivided 
into 2 and 3 sub-clumps respectively. Collectively, these 
sub-clumps had properties very similar to their respec- 
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Table 3. Properties of the 3D-clumps 
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Fig. 5. SCUBA map of the 850 fim continuum emission 
(colour scale) oveiplotted with all the J=l-0 3D-clumps 
integrated intensity maps (contours and numbers). The num- 
bering of the cores is based on their peak intensity. The differ- 
ent contour style and colours identify clumps starting at dif- 
ferent contour levels. The white solid contours are the clumps 
with the stronger integrated emission, with contours starting 
at 1 K kms"' , the dashed yellow contours start at 0.3 K kms"' , 
and the solid red contours start at 0.1 Kkms ' step. For all 
clumps, the contour step is 0.1 Kkms"'. 



live C'^O clumps. The presence of these possible sub- 
clumps does not significantly alter the interpretation of 
the region for the purpose of our analysis, indicating 
that the clumps adequately describe Serpens. 



3.3. Structure of the region: combining information 
from 2D and 3D-ciumps 

The north region has two clear clumps unequivocally 
identified in both 2D and 3D methods: 2D-clumps 
A and B, which correspond to 3D-clumps 1 and 2 
respectively. Both peak close to the position of the 
strongest submillimetre sources in this region (SMMl 
and SMM9), and trace the gas around them in good 
agreement to the cold dense dust traced by the 850 i-im 
emission. 

A region with higher velocity gas was detected 
with the 2D analysis as 2D-clump C which corre- 
sponds to 3D-clump 5. This region has quite strong in- 
tegrated emission making it detectable in the 2D search. 
However, as it peaks at a very similar velocity to clump 
1/A, the 3D search failed to separate these two in some 
of our trial runs of the 3D analysis. This clump is as- 
sociated with very little submillimetre continuum emis- 
sion but quite strong C'^O (and C"^0) emission. The 
fact that it is also seen in N2H+ (Olmi & Testi 2002) and 
not in '^CO tracing outflows (SFLPaper), is consistent 
with the possibility of this being a denser region, close 
to being bound, directly associated with the NW sub- 
cluster. It could, for example, be a very young prestel- 
lar core about to become gravitationally unstable and 
collapse (Walsh et al. 2007). 

A region detected with the 3D analysis which was 
not seen in the 2D search was 3D-clump 8. This clump 
is detected at high velocities (8.8 kms"') and seems to 
surround the clump 2/B associated with SMM9, per- 
haps as a shell. Although apparently somewhat super- 
virial clump, if aff'ected by a depletion of C"^0 by a 
factor of 2-3 (Appendix B), this clump could be gas 
undergoing gravitational collapse. 

Finally there is also a low velocity region situated 
at the left of the main clumps of the NW sub-cluster - 
detected as a single clump with the 2D method (clump 
D, on Fig. 4) and as three separate clumps with the 3D 
CLUMPFiND (11, 13 and 15 on Figure 5). This region has a 
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very small mass (~ 1 - 2 M©) and is about 5 times super- 
virial. It seems to be a quiescent region at lower veloc- 
ities than the main cloud and connecting to the main 
cloud very close to the edge of the NW sub-cluster as 
seen on dust emission. 

The bulk of emission on this NW sub-cluster 
presents a very coherent structure in space and veloc- 
ity throughout. It does not appear to be as filamentary 
as the SE sub-cluster and the emission appears confined 
to relatively dense, cool compact regions. 

The SE sub-cluster is quite different from the NW 
cluster, both in spatial structure and velocity, even 
though this is not obvious from the dust emission. Note 
the higher virial ratios for the clumps within the SE sub- 
cluster, when compared to the ones in the NW, support- 
ing once again the idea of more kinetic support in the 
south, even when the 3D velocity- separated clumps are 
considered. Comparing the 2D and 3D results, shows 
the emission is more complex with none of the 

gas emission peaks coincide with any of the compact 
submillimetre sources. The main peaks of the C'^O 
emission in this region lie in the filament seen in dust 
continuum emission, between the compact sources. The 
2D-clumps E and F were detected as 3D-clumps 10 and 
3 respectively. However, the 3D search found a more 
difliise clump, 3D-clump 12, which peaks east of the fil- 
ament but with its edges still overlapping spatially with 
3D-clump 10 and 3, having lower velocities than these 
two: 7.36 kms"' of clump 12, versus 7.68 kms"' and 
8.67 kms"' of clump 10 and 3 respectively. Despite be- 
ing adjacent, and with overlapping edges, 3D-clumps 
3 and 10 have a difiference of 1 kms ' between their 
peak velocities. There is a similar velocity difference 
between clump 10 and clump 6, west of the filament: 
3D-clump 6, detected in the 2D analysis as 2D-clump 
G, hasapeakvelocity of 8.51 kms"^ ~0.8kms"' higher 
than its neighbour. These four 3D-clumps (3, 6, 10 
and 12), with two sets of diff'erent peak velocities (at 
~7.5 kms"^ and ~8.5 kms"'), overlap with each other 
at low intensities mainly throughout this filamentary 
structure of the SE sub-cluster, even though their emis- 
sion peaks are spatially offset. This also shows that the 
double velocity structure in the SE sub-cluster (Section 
2.1.2) is to some extent, recoverable from a single line 
fit using a 3D clumpfind analysis. 

The remaining clumps detected in the SE sub- 
clusters trace the less dense gas around this main fila- 
ment. These were not detected in the 2D search mainly 
due to their very narrow line widths, between 0.3 and 
0.5 kms"', making them faint in integrated intensity 
maps. Note that the dominant emission detected east 
of the filament has lower velocities (3D-clump 7 has 
a peak velocity of 8.09 kms"'), whereas the regions de- 
tected to the west have higher velocities (3D-clumps 
4, 9 and 16), with mean velocities from 8.20 kms~^ to 
8.80 kms ^ 

Globally, there appears to be a velocity gradient 
from east to west of nearly 1 kms~' over slightly more 
than 0.1 pc. However, this is not a smooth gradient 
throughout, as in the filamentary structure there are 
spatially-overlapping clumps with very different veloci- 



ties. This velocity structure is further investigated using 
the C^^O lines, which are not split by hyperfine struc- 
ture in Section 4.2. 

4. Dynamics of the cloud: velocity and line 
width 

4.1. Outflows 

One important issue when studying gas dynamics in re- 
gions of active star formation is the extent to which the 
fine widths of molecular species are influenced by out- 
flows. Using the available data, we looked for the in- 
fluence of outflows on the size scale of the cores by 
investigating the spectra associated with all the submil- 
limetre sources, looking for possible wing emission. 

Although wings on C^^O lines have proven to be 
able to trace outflow interaction (Fuller & Ladd 2002), 
in Serpens and with the 0.45K rms noise of our dataset 
(Section 2.1) no wings were found. The lines towards 
sources with known outflows are well fitted by a sin- 
gle Gaussian. For example. Fig 3 (top panel) shows the 
C'^O towards SMMl, a source known to have an out- 
flow, e.g. Hurt & Barsony 1996, is well represented by a 
single Gaussian component. In the SFLPaper, we have 
also searched for evidence of the influence of outflows 
in the C^^O emission by comparing the C'^O J=3— »2 
emission to the '^CO J=3^2 emission tracing the out- 
flows. No correlation nor anti-correlation between the 
C'**0 emission and the outflows is found in the re- 
gion. Both these approaches lead us to conclude that 
the Serpens C^^O emission is not influenced by out- 
flows and, therefore, the velocity components we detect 
in C^^O are related to the global cloud dynamics. 

4.2. Position - Velocity structure 

As revealed by flie C^^O (Section 3.1), the NW region 
is mostly traced by higher velocity emission, the ex- 
ception being the region offset to the east of the sub- 
cluster. On the other hand, the SE is not so homoge- 
neous, containing both higher and lower velocity com- 
ponents, which overlap approximately where the fila- 
mentary structure is seen in the continuum observa- 
tions. The C^^O shows very well defined double peaked 
emission along the southern region, that starts disap- 
pearing as we move north. Figure 3 shows two exam- 
ples of spectra in this region, and Fig. 9 shows the evo- 
lution of the double component along the map. The ex- 
istence of a double-peaked spectrum in other optically 
thin tracers such as N2H^ (Olmi & Testi 2002) rules out 
self absorption as an explanation of the double peaked 
C'^O emission (cf Section 2.1.2). 

On the basis of a study of the fine centroid veloc- 
ity (despite the presence of double peaked lines) Olmi 
& Testi (2002) argued that the region is undergoing 
global rotation. However position- velocity diagrams of 
horizontal slices along the C'^O map are incompati- 
ble with this interpretation (Fig. 6 and Fig. 7 to 9). 
Moving from north to south, and slicing at the declina- 
tion of each SMM source, we can see the two separate 
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1*1 6'M' 



60CO 



PV1 
PV2 



i1*14'Cn" 



1*12'CJ' 




Fig. 6. SCUBA 850 map of Serpens in colour scale and 
contours (contours as in Fig.l), showing the cuts used for the 
position-velocity diagrams of C'*0 J=l-0 emission (Fig. 7, 8 
and 9). The cuts made crossing the SMM sources are denoted 
as PV#, where # is the number of the cut, starting from the 
north. There are two cuts which are denoted as "dust" which 
do not intercept any source, but were made to understand the 
velocity structure around the dust filament. 



clouds, very well distinguished close to SMM 11. For 
simple rotation, we would expect to see a smooth gra- 
dient along the velocity axis as the RA changes. Instead 
we observe two velocity components, clearly separated 
in the southern part of Serpens (see e.g. PVIO) and 
merging together when moving to the north of the sub- 
cluster (see e.g. PV7). At this point, the two compo- 
nents are barely distinct lines, producing broad, non- 
gaussian profile. Furthermore, the SMM sources in the 
SB sub-cluster appear at the edges of the double veloc- 
ity region (hereafter referred to as the interface), whilst 
the filamentary structure seen in dust follows the inter- 
face region itself (see the PV diagrams labeled as "dust" 
on Fig. 8 and 9). 

The lower velocity cloud (hereafter LVC) appears 
to be interacting with the high velocity cloud (here- 
after HVC), apparently provoking the enhanced dust 
emission between SMM2, SMM3, SMM4 and SMM6 - 
and also the elongated filament that extends south to- 
wards SMM 11 and beyond. A dynamical interaction 
between two clouds, as indicated by this space-velocity 
structure and the turbulent motions found towards the 
south, might have triggered this episode of star forma- 
tion along the filament (Sec. 6). 



PVl 




PV2 




PV3 




py4 ,„.wi,,. =n„i 

Fig. 7. Position-velocity diagrams of the C'*0 J=l— *0 emis- 
sion. The cuts are horizontal slices of the map, as shown 
in Fig. 6. The sources name and RA position are indicated 
in each figure (dashed line). The declinations are presented 
in Table 1 and the RA varies from 18''30"'06' to 18''29"'46' 
(from to -300 " offset respectively). The gray scale shows 
the line intensity in T^. 
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Fig. 8. Same type of position- velocity diagrams as Fig.7 for Fig. 9. Fig. 8 continued: remaining position- velocity dia- 

the SE sub-cluster. The RA also varies from IS'-SO'^Oe' to grams of the SE sub-cluster as plotted on Fig. 6. 
18''29""46' (from to -300 " offset respectively). PV diagrams 
displayed in descending declination, as they appear on Fig. 6. 

The data were first rebinned to 0.1 kms"' veloc- 

4.3. Decomposition of the C^'O line components ity channels. Then, for each spectrum the line was fit- 

ted with a single Gaussian and then a double Gaussian. 

To investigate the velocity structure of the C'**0 J=l— >0 The two Gaussian fit was selected as the model for the 

emission we have decomposed the datacube in two, by line only if: i) the difference between the central ve- 

fitting two velocity components to the C'^O spectra and locities of the two Gaussian fit (AV) was greater than 

then creating a model datacube from the Gaussian fits, 0.35 kms"' or ii) both lines were relatively strong with 

for each of the two components. The results of this de- the peak intensity ratio of the stronger to the weaker 

composition are shown in Fig. 10 and Fig. 11. fine less than 2.4. The value of 2.4 was determined by 
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Fig. 10. Integrated intensity maps as a result of the separation of tiie two line components of the J=l-0 transition, under 
the assumption of two different clouds seen along the line of sight. The LVC is shown on the left and HVC on the right. The 
background grey scale shows the 850 yum emission tracing the cold dust, with its respective submillimetre sources plotted as 
triangles. The contours represent the integrated intensity of the modelled Gaussians (fitting the data in T^). Contour key: left 
(LVC) at 0.2 Kkms-' (dashed) and 1.5, 2.5, 3.2 and 4.0 Kkms"' (solid); right (HVC) at 2.5 and 3.0 Kkms"' (dashed) and 3.2, 
4.0, 5.0, 5.5 and 6.0 Kkms"' (solid). 




fit (. Ht 

Fig. 11. Velocity structure maps of Serpens, as a result of the separation of the two line components of the C"*0 J=l-0 
transition. As in Fig. 10, the LVC is represented on the left and HVC on the right. The submillimetre sources are plotted as 
triangles, and the contours represent the integrated intensity as in Fig. 10. The colour scale is now the centroid velocity of the 
same modelled Gaussians, where the light pink colour represents the lack of a fit to that velocity component. 



a careful analysis of various line fits which showed that 
if the ratio was more than 2.4, the weaker line fit was 
poorly fit. The remaining spectra were fitted with a sin- 
gle Gaussian. An example of this fitting is shown on 
Figure 3 (lower panel). From this fitting procedure, two 
model datacubes were created, one for each velocity 
component, allowing us to study the two clouds sep- 
arately. The higher velocity component (HVC) of the 
double peaked lines, as well as the single lines with 
central velocity greater than 7.8 kms"', were included 
in the HVC datacube; lower-velocity lines and single 
lines peaking below 7.8 kms ' were incorporated in the 
LVC model datacube. 

Figure 10 shows the spatial distribution of the LVC 
and HVC using the integrated intensity from the model 
datacubes. Overall, the HVC traces the distribution of 
the 850 jum continuum emission better than the LVC. 



The HVC emission is stronger in the north, but it lies 
along the filament containing both sub-clusters, extend- 
ing roughly in a SE-NW direction. The LVC is roughly 
aligned along the S-N direction and is stronger in the 
south, where it meets the HVC. 

Figure 11 shows again the integrated intensity 
of the modelled datacubes for each component, but 
shows also the velocity structure of each cloud. In the 
NW sub-cluster, the HVC appears at velocities around 
8.4 kms"', with the exception of a few regions at the 
edges of the cloud, which appear to reach velocities as 
high as 8.8 kms"'. The region which stands out from 
the bulk of this sub-cluster is the region SW of SMMl, 
which is not present in the dust emission even though 
it is rather strong in gas emission, having the highest 
velocities of the entire cloud (reaching 9 kms"'). The 
LVC, in the NW sub-cluster, is spatially offset to east. 
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with velocities of 7.5 - 7.8 kms"', similar to most of the 
emission in the south. 

The region between the two sub-clusters, dominated 
by the emission from the HVC, has the systemic veloc- 
ity of Serpens (around 8.0 kms"') possibly due to the 
merging of the two components. Note that the emis- 
sion here is rather weak, and the presence of SVS2, a 
more evolved (flat spectrum) near-IR source (Kaas et al. 
2004), suggests this region may be more evolved. 

In the SE sub-cluster, the HVC velocities range 
from 8 and 8.5 kms"', being higher towards the south- 
ern end of the filament. On the other hand, the LVC 
shows a velocity gradient increasing from west to east - 
contrary to the HVC. The material west of the south- 
ern filament, has velocities of about 6.8 - 7 kms"'. 
At the centre of the filament the velocities are around 
7.5 kms"', translating into a gradient of ~ 5 kms^'pc"'. 
To the east of the filament the velocities are approx- 
imately constant and around 7.5 kms ^ Therefore, it 
seems that the clouds have a greater offset in velocities 
in the far-south end of the filament, converging into one 
intermediate velocity as one moves north. When two 
lines can no longer be separated, the emission becomes 
a single broader line, centered at the intermediate ve- 
locities (~ 8 kms"'). 

The line width in the SE sub-cluster, specially 
where the two components merge, is around 2 kms"'. 
This is almost twice the line width of the NW 
(~ 1 kms"'). This difference is reflected as a four times 
higher kinetic support in the SE region, in compari- 
son to the NW region. This is consistent with the C'^O 
J=l— >0 analysis (Section 3. 1), where it was showed that 
the NW sub-cluster is a bound structure, whereas the 
SE sub-cluster was somewhat super-virial. The SE is 
therefore much more dynamic than the NW, as already 
foreseen by the C'^O analysis. 



5. Physical properties: Temperatures and 
column densities 

The optical depth of the C"*0 can be estimated from the 
ratio of the integrated intensities of the C'^O and C'^O 
J=l— »0 transitions. Over the mapped region the ob- 
served ratio shows little coherent spatial structure and 
is approximately constant with a value consistent with 
the abundance ratio of the species, ~3.5 (e.g. Penzias 
1980; Frerking et al. 1982) implying the emission from 
both species is optically thin. 

To better understand the correlation between the 
dust and gas in this region. Figure 12 shows a pixel- 
by-pixel comparison of the 850 //m flux density against 
the integrated intensity of the three transitions of C"*0, 
all convolved to a common resolution of 24". For the 
purpose of these scatter plots, we have oversampled the 
data to a pixel size of 2.5", in order to better distinguish 
the trends. 

Overall, the distribution of points is very similar for 
the three transitions. There is a general correlation be- 
tween dust and gas, especially for the weaker emission 
(Fig. 12). However, the distributions also show structure 
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Fig. 12. Scatter plots of the SCUBA 850 fim dust emission 
against the C"*0 integrated intensity over the whole range 
of velocities (top: J=l— >0, middle: J=2— >1, bottom: J=3— *2). 
Blue points show the NW sub-cluster and green points the 
SE sub-cluster. The four positions chosen in each sub-cluster 
to investigate with non-LTE modelling shown on the middle 
panel for the NW region and the lower panel for the SE re- 
gion. These positions are also indicated on Fig. 13. 
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Fig. 13. Map showing the positions of the selected regions 
for the non-LTE RADEX study indicated by blue and green 
circles (for the positions in the NW and SE sub-clusters, re- 
spectively) and labeled as in Fig. 12. The contours and colour 
scale show the SCUBA 850 /jm emission as in Fig. 6. 



which consistently appears across all three transitions. 
The very prominent peaks of dust emission correspond- 
ing to the stronger submillimetre sources are obvious, 
and although in general there is an increase in the C'^O 
emission at these positions, the dust peaks do not cor- 
respond to global peaks in the emission. Indeed 
the nature of the relationship between the C'**0 emis- 
sion and the dust appears different in the NW and SE 
sub-clusters. 

Focusing on the NW region (blue in Fig. 12), the 
plots are dominated by the two dust peaks, each of 
which is associated with a well defined, but separate, in- 
crease in emission. Comparing the C'^O intensity, 
the emission becomes weaker moving to higher energy 
transitions. On the other hand, the SE sub-cluster (green 
in the figure) shows a different trend from transition to 
transition, becoming stronger at higher transitions. In 
addition, there appears to be a more pronounced gen- 
eral correlation in this region between the dust and line 
emission. Neveitheless, there are clearly structures de- 
parting from this trend: several 850 fj.m peaks corre- 
sponding to SMM sources; and C'^O peaks, which do 
not have significant submillimetre emission. 

5.1. LTE Analysis 

From the dust continuum emission, the volume den- 
sities in the Serpens sub-clusters are typically higher 
than the critical densities of each of the three transitions 
(Table 4) observed here. We therefore initially calcu- 
late the gas properties assuming LTE (local thermody- 
namic equilibrium) using a rotation diagram analysis. 
Despite its uncertainties, the rotation diagram method 
is robust in retrieving the column density structure and 
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Fig. 14. LTE excitation temperature map in colour scale. The 
dotted black contours show the dust 850 yum emission at 0.6, 
1.2 and 1.8 Jy beam"'. 
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Fig. 15. Column density map (colour scale) derived from 
the rotation diagram method. The dust 850 fim emission is 
overplotted in dotted black contours, at the same levels as in 
Fig. 14. 



trends throughout the region, as well as the approximate 
absolute column densities. 

Fig. 14 shows a map of the excitation temperature 
across the region constructed from the 24" resolution 
integrated intensity maps sampled with 5" pixels (as 
in the original IRAM data). This shows the NW and 
SE sub-clusters to have different temperature structures. 
The NW appears very homogeneous with no significant 
temperature peaks and with temperatures ranging from 
9 to 10 K. In contrast the SE region has both higher 
temperatures, ranging from ~ 10 to 14 K and a much 
more peaked distribution. Interestingly, this enhanced 
temperature in the south does not peak on the SMM 
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Table 4. Critical densities (ticnticai) at lOK and 20K for C'^O 





Au 


Ku 


Identical 


transition 


(s->) 


(10-" cm^s-i) 


(10^ 






lOK - 20K 


lOK - 20K 


J=1^0 


6.266x10-** 


3.3-3.3 


1.9-1.9 


J=2^1 


6.011x10"^ 


7.2 - 6.5 


8.3-9.3 


J=3^2 


2.172x10-*' 


7.9-7.1 


27-30 



These values were derived using the information provided by the Leiden Atomic and Molecular Database^ 
(Schoeier et al. 2005): Au being the Einstein Coefficient for the upper level, and Ku the respective collision rate at 
both 10 K and 20 K. 



protostars but rather between them, along the dust fila- 
ment which corresponds to the interface region seen on 
the PV diagrams (Fig. 8 and 9). 

The C'^O column density map (Fig. 15) calculated 
from the rotation diagram recovers more of the dust 
structure than the temperature map. Both the south and 
north sub-clusters are evident as denser regions, even 
though the dust and gas column densities peaks are 
not always coincident, especially in the SE. The mean 
C'^O column densities are very similar in the north and 
the south. The regions with higher gas column density 
(the entire NW sub-cluster and the filament between 
SMMl 1 and SMM2 in the SE sub-cluster) have a lower 
temperature. Conversely the regions with slightly lower 
gas column density (between SMM2, SMM4, SMM3 
and SMM6 in the SE sub-cluster) have higher tem- 
perature. The region south-west of the NW sub-cluster 
which appears to have a relatively high gas column den- 
sity seems to have very similar properties to the rest of 
the NW sub-cluster and yet it is not detected in dust 
emission. 

5.2. Non-LTE Analysis 

To better understand the physical conditions of Serpens 
and the apparent discrepancies between the gas and dust 
emission, we have selected 4 key positions in the NW 
and another 4 in the SE using the scatter plots (Fig. 12) 
for more detailed analysis. The selected positions are 
indicated on Fig. 13. These positions correspond to in- 
teresting features in the correlations between the dust 
and gas emission, selected to span the range of the cor- 
relation. For these positions we performed a non-LTE 
analysis using RADEX"'. 

We created a 500 x 500 grid of gas column den- 
sity (ranging from 10'^ to 10'^ cm"^) and temperature 
(ranging from 5 K to 40 K). For each grid point we 
used RADEX to calculate the C^^O integrated intensi- 
ties for all three transitions (denoted as For 
each position modelled, a volume density determined 



' http://www.strw.leidenuniv.nl/moldata/ 
RADEX is a statistical equilibrium radiative transfer 
code, available as part of the Leiden Atomic and Molecular 
Database. The formalism adopted in RADEX is summarized 
in van der Tak et al. (2007) 



from the submillimetre dust continuum emission was 
adopted, assuming a cloud depth of 0.2 pc (based on 
the projected size of the dust emission), a dust opac- 
ity of 0.02 cm^g-i at 850 jum (van der Tak et al. 1999; 
Johnstone & Bally 2006) and a dust temperature of 
10 K for all but three positions. The three exceptions 
are: position NA (=SMM1) where 38 K was adopted 
from the SED fit by Davis et al. (1999), and positions 
NB(=SMM9) and SA(=SMM4), where we adopted a 
temperature of 25 K, consistent with the > 20 K deter- 
mined by Davis et al. (1999). For each of the southern 
positions, we assumed the H2 volume density to be the 
same for both C^^O velocity components. Changing the 
dust temperature or the assumed cloud depth changes 
the estimated the H2 volume densities, however this 
only becomes important if the derived volume densi- 
ties becomes lower than the critical densities for our 
transitions. We have tested these effects using RADEX, 
and the resulting column densities and kinetic tempera- 
tures remain unaffected by changes in the assumed dust 
temperature between 10 K and 40 K, or in the assumed 
depth between O.lpc and 0.3pc. If the ttansitions are 
thermaUsed, only the fractional abimdance of C'^O will 
be affected by changing the assumed H2 column den- 
sity. 

The central velocity, Une widths and integrated in- 
tensity of each transition were retrieved from the data 
by fitting the average spectrum within a 5" radius of 
each position. The central velocity and line widths 
shown in Table 5 are the average over the three transi- 
tions, and have an uncertainty of the order of 0. 1 kms'^ 
For the 4 positions in the NW sub-cluster, this pro- 
cedure is straight-forward as the lines of all three 
ttansitions are well represented by single Gaussians. 
However, the spectra of the SE sub-cluster positions, 
having two velocity components, was separately fitted 
with 2 Gaussians in order to investigate any possible 
differences between the two components. We used z.x^ 
comparison to find the best fit of the RADEX models to 
the observed ratios /(i_o)//(2-i) and /(i_o)//(3-2) as well 
as the absolute value of /(i-o). 

In calculating the the relative errors of the in- 
put integrated intensities were assumed to be equal to 
the relative errors in the intensity, Tmb), given by the 
rm.s. of the line fit given by class. The results from the 
RADEX models are presented in Table 5. For compari- 
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Table 5. Modelling results for the 8 positions selected from the scatter plots (App. A). 





TT 

H2 


Line 


Central 


non-LTE 


non-LTE 


T 'I'L' 


• * 

Position 


volume density 


width 


velocity 


column density 


tkin 


T 

exc 




(x 10^ cm"^) 


(kms ) 


(kms ) 


(x 10 cm ) 


(K) 


(K) 


NA 


2.60 


1.6 


8.5 


15.4 


11.7 


10.3 


NB 


2.25 


1.5 


8.4 


11.9 


10.6 


9.5 


NC 


1.62 


1.9 


8.5 


11.9 


10.5 


9.5 


ND 


6.50 


1.4 


8.5 


13.5 


11.4 


9.5 


SAl 


2.98 


2.2 


7.8 


6.4 


18.2 


13.4 


SA2 


2.98 


1.0 


8.3 


4.2 


6.6 




SBl 


3.99 


1.5 


6.9 


6.0 


11.9 


12.7 


SB2 


3.99 


1.1 


8.5 


4.4 


14.8 




SCI 


6.28 


1.7 


7.7 


5.8 


11.7 


10.2 


SC2 


6.28 


1.2 


8.7 


1.4 


14.2 




SDl 


7.81 


1.1 


6.9 


7.3 


16.6 


12.8 


SD2 


7.81 


1.4 


8.2 


5.8 


14.4 





Positions starting with N are in the NW sub-cluster while positions in the southern region start with S. For the 
positions in the SE, the labels 1 and 2 identify the low and high velocity components respectively. 



son, the table also includes the temperature and column 
densities retrieved using the rotation diagram method. 
The surfaces are shown in App. A, while App. B 
shows the input and best fit integrated intensities for the 
three Unes, as well as the impUed abundances from the 
non-LTE analysis. 

Table 5 shows that in the NW region the LTE 
(rotation diagram) and non-LTE (RADEX) analyses 
are in good agreement. Both the rotation diagram and 
RADEX show variations in Tex of only 1 K. They pro- 
duce absolute values of temperature which differ by at 
most 10%, and the trend in temperature between posi- 
tions is similar. Comparison of the C'^O column den- 
sities and the H2 column density calculated from the 
dust continuum emission at the same positions, imply 
C'^O abundances a factor of ~2.5 smaller than typi- 
cal values (App. B). Given the low optical depth of 
the C'^O emission (Sec. 5), the low C'^O column den- 
sity, and hence abundance, implies that the C'^O is de- 
pleted with respect to the molecular hydrogen, even in 
the warmer envelope of SMMl, consistent with the re- 
sults of Hogerheijde et al. (1999). 

Somewhat surprisingly given the low optical depth 
of the C'**0 and despite its depletion, neither the LTE 
nor the non-LTE analysis finds evidence of increased 
temperatures towards the apparently warmer inner re- 
gions of the embedded protostars. Although the dust 
emission indicates the presence of warm dust towards 
the protostars (Davis et al. 1999), the C"^0 emission 
implies low and uniform temperatures. Although di- 
lution of the warm inner region within the 24" beam 
may contribute to the difficulty in detecting the warmest 
gas, it is surprising that no evidence of any tempera- 
ture increase is seen. CO is predicted to freeze-out on 
to grain surfaces at temperatures below ~ 18K, consis- 
tent with the low C^^O excitation temperature, but not at 
dust temperatures of, for instance, ~ 30K seen towards 
SMMl. Since this is above the sublimation temperature 



of pure CO ice it is possible the CO could be trapped in 
a water rich ice, which would only sublime and return 
CO to the gas phase at temperatures of ~100 K (e.g. 
Visser et al. 2009). 

The southern region is more complex. In terms of 
column density, for all four positions the non-LTE re- 
sults show the LVC to have sUghtly higher column den- 
sity that the HVC. We find that the LTE and non-LTE 
approaches agree in the sense that the northern posi- 
tions and SD have higher values of C'^O total column 
density (summed over both velocity components where 
necessary). These are followed with decreasing column 
density by SB, then SA and finally SC . 

SA is at the position of SMM4, where there are two 
components of the C^**0 emission, a strong low veloc- 
ity component (SAl) plus a weaker high velocity com- 
ponent (SA2). In the J=3— >2 transition, the high ve- 
locity component becomes faint and difficult to sepa- 
rate from the lower velocity component (see Table B.l 
in App. B). This weak J=3— >2 emission constrains the 
temperature to 6.6 K for the higher velocity compo- 
nent (SA2). At SC the two cloud components are also 
significantly blended. The temperature of the LVC at 
this position (SCI) is ~11 K, constrained within ~1 K 
(App. A), with the weaker HVC warmer, but somewhat 
less well constrained. 

In general, in the south, the lower velocity compo- 
nent has a higher temperature toward the most central 
positions studied (SA and SD). Then at SB and SC, 
at the edges of the dust emission, the temperatures are 
similar with both components still higher than the tem- 
peratures generally found to the north (~12 K and 14 K 
in the south versus ~11 K in the north). Overall, and 
with the exception of SA, the HVC has higher temper- 
atures than in the north, around ~14 K. On the other 
hand, the LVC traces the temperature trend as identi- 
fied by the LTE study (Fig. 14) better than the HVC, 
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but the absolute LTE temperatures are between the non- 
LTE values for LVC and HVC. 

Therefore, we conclude the temperature rise toward 
the south is real. Such a rise is consistent with a scenario 
where this region is tracing the interaction/collision be- 
tween two clouds, with a shock layer with higher tem- 
peratures and complex motions at the interface. 

6. Discussion 

6.1. Two different sub-clusters in Serpens 

Our study of the Serpens Main Cluster has shown that 
two apparently very similar protoclusters as seen in sub- 
millimetre dust continuum emission can reveal very dif- 
ferent dynamical and physical properties in molecular 
lines. Despite aU the outflows seen in '^CO in the re- 
gion, the denser gas around the cores seen in C^^O and 
C'^O does not seem to be perturbed and is able to pro- 
vide details of the quiescent material in the cloud. 

In the NW sub-cluster the bulk of emission has a 
velocity around 8.5 kms"'. However, there is a lower 
velocity component of the gas east of the sub-cluster 
(Fig. 7) with the transition between these component 
being rather smooth. The velocity difference between 
the submillimetre sources in this sub-cluster is small, 
ranging from 0.1 to 0.3 kms"'. 

The physical conditions in this NW sub-cluster are 
also rather coherent. Temperatures and column den- 
sities derived from both LTE and non-LTE analyses 
are consistent and show little variation within the sub- 
cluster. The C'^^O emission peaks are mostly consistent 
with the dust peaks. Clump-finding studies of this re- 
gion retrieved two main peaks which are directly re- 
lated to the two stronger submillimetre sources in the 
NW sub-cluster: SMMl and SMM9. However there are 
no evident temperature peaks associated with the sub- 
millimetre sources. The remaining gas emission in the 
sub-cluster is either associated with these main peaks or 
weaker structures surrounding the main bulk of the dust 
emission. The gas column density very closely follows 
the clumps/integrated intensity distribution of the gas, 
particularly in the lower J transitions tracing the colder 
gas. 

The SE sub-cluster on the other hand is a much 
richer region in its dynamics and properties. There 
are two velocity components/clouds along the line of 
sight, clearly identified using both clump-finding and 
position-velocity diagrams. These two clouds appear 
to be interacting. They are more ofi^set in velocity in 
the south and start to mix moving to the north within 
the sub-cluster. Furthermore, the submillimetre sources 
in this region appear at the edges of the interface of 
the two components, whereas the dust filament ap- 
pears in the interface. Most of the southern submillime- 
tre sources appear to have a stronger association with 
the HVC, despite having some emission from the LVC 
along the same line of sight. A counterexample how- 
ever is SMM2 which, as can be seen in the PV dia- 
grams, has a stronger C'**0 1=1— >0 lower velocity com- 
ponent. The overall dust filament, as seen in 850 fim. 



coincides with the N-S lane where the two components 
overlap, suggesting it is tracing the interface region be- 
tween the components, the region where they are in- 
teracting. Ultimately, this interaction might have been 
responsible for triggering the star formation episode in 
the SE sub-cluster. 

In contrast to the NW sub-cluster, the LTE tem- 
perature in the SE sub-cluster is both higher and more 
structured, peaking close to the ridge of dust contin- 
uum emission. Unhke the north, the two velocity com- 
ponents in the south are difficult to fit with a single well 
defined temperature. The general trend, however, points 
to higher temperatures in the southern sub-cluster than 
in the northern sub-cluster. 

The modelled column density map (Fig. 15) closely 
traces the emission from the lower transitions (J=1^0). 
The high C'^O column density regions in the SE are 
not associated with any of the submillimetre sources, 
but rather the southern filament. The region with en- 
hanced temperature, however, does not coincide with 
the highest column density regions. The uniform dust 
emission over this SE region results from the southern 
filament having lower temperature but higher column 
density whereas the northern part of the SE sub-cluster 
is slightly less dense, but warmer, resulting in equiva- 
lent 850 jum dust emission. 

6.2. A Proposed Scenario 

The velocity, temperature and density structure of 
Serpens suggest a more complex picture than simple 
rotation which has previously been invoked to explain 
the velocity structure (e.g. Olmi & Testi 2002). 

It is known that cloud-cloud or flow collisions hap- 
pen in the Galaxy as molecular clouds move within 
the spiral arms. Furthermore, simulations of cloud- 
cloud collisions (e.g. Kitsionas & Whitworth 2007) 
have shown that density enhancements in the collision 
layers can be high enough to trigger star formation. 
Additionally, clouds are commonly seen as filamentary 
structures, not only during, but also prior, to star for- 
mation. We suggest that the two velocity components 
seen in Serpens are tracing two clouds along the line of 
sight and that the interaction of these clouds is a key 
ingredient in the star formation in Serpens. 

We propose that we are seeing two somewhat fila- 
mentary clouds traveling toward each other and collid- 
ing where the southern sub-cluster is being formed. The 
cloud coming toward us is to the east while the cloud 
moving away from us is to the west, and represents the 
main cloud. An inclination angle between the two fila- 
ments could explain why the two velocity components 
are spatially ofl'set in the north but overlapping in the 
south. This scenario explains both the double peaked 
profiles of the optically thin fines and their distribution 
along what has previously been identified as the 'rota- 
tion axis' of this region. 

If the north region was initially close to collapse, 
the direct collision of the clouds in the south could in- 
directly trigger or speed up this collapse in the north 
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without significantly enhancing the temperature or per- 
turbing the intrinsic, 'well behaved' velocity and col- 
umn density structure. In the south, however, such a 
collision makes it easy to understand why the density 
and temperature enhancements are not necessarily as- 
sociated with the sources, as they are being generated 
by an external trigger: the colhsion. 

Note in addition, that in the south, unlike the major- 
ity of the sources in the north, there is a poor correlation 
between the submillimetre sources (Davis et al. 1999) 
and 24 /zm sources (Harvey et al. 2007b), as shown in 
Fig. 1, suggesting a wider spread of ages of the proto- 
stars in the south than in the north. Such an age spread 
would be consistent with a collision in the sense that 
a collision is not a one-off event but rather an ongoing 
process. 

A first test to this collision scenario is the timescale 
for which such clouds would cross each other, their in- 
teraction time. We assume that each cloud is a filament 
of radius of 0. 1 pc (similar to the size of the dust 850 fim 
emission). In addition we adopt a collision velocity of 
I km s ' (approximately the mean observed velocity 
difference, along the line of sight, between the two com- 
ponents). The timescale from when the clouds start col- 
liding until they are completely separated, assuming a 
head on collision, is 4 x 10^ years, consistent with the 
estimated ~ 10^ year age of the region (Harvey et al. 
2007a; Kaas et al. 2004). 

Several simulations of cloud collisions such as pro- 
posed here exist in the hterature. For example, SPH 
simulations of clump-clump collisions from Kitsionas 
& Whitworth (2007), have shown that two approach- 
ing clumps with a slow collision velocity of 1 kms ' 
(Mach number of 5), can indeed trigger star formation 
in the colhsion layer. Specific simulations of the pro- 
posed collision in Serpens will be presented in a subse- 
quent paper (Duarte-Cabral et al. 2010 in prep.) 

7. Conclusions 

The dynamics and structure of the Serpens Main 
Cluster have been studied in detail, as an example of 
a complex low mass star cluster forming region. This 
study has provided a view of the dynamics and structure 
of the region. The Main Cluster comprises a very young 
star forming region subdivided into two sub-clusters. A 
careful investigation of the clump structure and excita- 
tion in the region shows that these two sub-clusters have 
similar overall masses but quite different properties. 

The NW sub-cluster is homogeneous in veloc- 
ity and temperature structure, with the submillimetre 
sources well correlated with the gas peaks as well as 
with the 24 ^um sources. On the other hand, in the SE 
sub-cluster there are two velocity components in the 
gas and the gas temperature is more variable. The gas 
column density and temperature peaks do not coincide 
with the submillimetre sources, but rather lay in the re- 
gions between them. Furthermore, the 24 /im sources in 
the south are poorly correlated with the dust emission. 

Our analysis suggests a scenario of cloud-cloud 
collision triggering the star formation in the SE cub- 



cluster, potentially inducing perturbations which indi- 
rectly affect the NW sub-cluster, hastening somewhat 
its collapse. SPH simulations of this scenario, to better 
understand how well it can reproduce observables, such 
as the velocity profile and the column density properties 
of the region, will be presented in a future paper. 
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Appendix A: surfaces for RADEX fits 

Figures A.l to A. 8) show the;^^^ surfaces from the non-LTE analysis of the Hne integrated intensity ratios for each 
of the positions studied. The^^ has been calculated using Eq. (A.l), where 1°^^ j is the observed integrated 

intensity of the transition between J„p and J/o^, A(x) is the uncertainty on the quantity x and, finally, /[y'^^y^ ) is 
the integrated intensity for each transition as modelled by RADEX, for each combination of temperature and gas 
column density. 



r/obs /robs _ mdex ; rradex \2 /robs /robs _ rradex / rradex 
'(1-0) '-'(2-1) ■'(l-0)'''(2-l) '(1-0)' -"(S-l) '(1-0)' '(3-2) 

. A(/(to)/^(ti)) J \ A^W^SP 



/ robs _ rradex 
'(1-0) '(1-0) 

A(/^-0)) 



(A.1) 



For each case, the;^'^ is plotted as a function of the RADEX output temperature and gas column density. Given 
the use of 3 quantities in the fit, we consider < 3, i.e. the reduced-;t'^ < 1 a good fit. All figures have contours 
dXx^ = 1, 2 and 3 with the exception of ND (Fig. A.4) and SD (Fig. A.8). 
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Fig. A.l. )^ surface for the integrated intensity ratios at position NA. 
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Fig. A.2.;jf^ surface for the integrated intensity ratios at position NB. 
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Fig. A.3.;^f^ surface for the integrated intensity ratios at position NC. 
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Fig. surface for the integrated intensity ratios at position ND. Contour at 25. 
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Fig. A.5. surfaces for the integrated intensity ratios at position SA: SAl (LVC) on the left, and SA2 (HVC) on 
the right. 
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Fig. A.6.;if^ surface for the integrated intensity ratios at position SB: SBl (LVC) on the left, and SB2 (HVC) on 
the right. 
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Col'jr-n dc'-'Mlv ' ' : Ctlorrn ccr'^ilv 

Fig. A.7.;tf^ surface for the integrated intensity ratios at position SC: SCI (LVC) on the left and SC2 (HVC) on 
the right. 
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Fig. A.8.;if^ surface for the integrated intensity ratio at position SD: SDl (LVC) on the left, and SD2 (HVC) on 
the right. Note the different colour scale for SC2. Contours are 6 for SDl and 18 for SD2 
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Appendix B: C^^O intensities and abundances 

Table B.l presents the best fit integrated intensities from the non-LTE (RADEX) modelling (§ A) , together with 
the observed values, and the implied C'^^O abundances. 

Table B.l. Modelled integrated intensities and resulting abundances. 







Observed 






Best fit 




H2 


c'»o 




Position 


1(1-0) 


1(2-1) 


1(3-2) 


1(1-0) 


1(2-1) 


1(3-2) 


column 


column 


fractional 






(Kkms-1) 






(Kkms-1) 




density 


density 


abundance 
















(1023cm-2) 


(10''*cm-2) 


(xlO-**) 


NA 


9.15 


10.65 


6.94 


9.23 


10.42 


7.28 


1.60 


15.4 


9.6 


NB 


7.47 


8.87 


4.76 


7.44 


8.15 


5.20 


1.39 


11.9 


8.6 


NC 


8.14 


10.35 


4.90 


8.13 


9.29 


5.57 


1.00 


11.9 


12.0 


ND 


7.90 


10.09 


4.82 


7.99 


8.88 


6.20 


4.01 


13.5 


3.4 


SAl 


4.69 


10.57 


8.46 


4.75 


9.84 


9.01 


1.84 


6.4 


5.8 


SA2 


2.42 


1.98 


0.50 


2.50 


2.02 


0.64 




4.2 




SBl 


4.63 


6.66 


4.20 


4.74 


6.72 


4.37 


2.46 


6.0 


4.2 


SB2 


3.32 


5.99 


4.31 


3.35 


5.73 


4.50 




4.4 




SCI 


4.61 


7.37 


3.87 


4.72 


6.82 


4.37 


3.87 


5.8 


1.9 


SC2 


1.17 


2.58 


1.45 


1.21 


2.23 


1.71 




1.4 




SDl 


4.81 


7.44 


7.31 


5.01 


8.53 


7.34 


4.82 


7.3 


2.7 


SD2 


4.41 


8.66 


4.66 


4.45 


7.42 


5.74 




5.8 





Positions of the NW sub-cluster are identified as starting with N, while the south positions start with S. For the positions in 
the SE, the labels 1 and 2 identify the lower (LVC) and higher velocity (HVC) components respectively. Where there are two 
velocity component lines, the C'^O fractional abundance was calculated using the total column density of C'*0, summing 
both components. 



For positions ND, SB, SC and SD, the H2 column densities derived from the dust and used to estimate the 
abundance of C"^0 were calculated using a dust temperature of lOK. Assuming a temperature of 15K for all 4 
positions (ND, SB, SC and SD) would reduce the H2 colunm densities by a factor of 2. 1, representing an equivalent 
rise of the fractal abundance of C^^O by the same amount. 

The derived C'^O fractional abundance (which is averaged along the line of sight) implies a depletion of 
C^^O of between a factor of 1.4 (for NC) and 4.3 (for SC), with an average of 2.5 compared to the abundance of 
1.7 X 10"^ in dark clouds (Frerking et al. 1982). Given that the ratio between C^^O and C"^0 has shown these 
two species to be optically thin, with an intensity ratio of ~3.5, a factor 2.5 depletion of C'^O implies the same 
depletion factor for C^^O. 
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